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, In the first part of these lectures we discuss the infrared limit of the spectrum of 

^ta* , the QCD Dirac operator. We discuss the global symmetries of the QCD partition 

Cn ' function and show that the Dirac spectrum near zero virtuality is determined by 

the pattern of spontaneous chiral symmetry breaking of a QCD-like partition func- 

Ch ' tion with additional bosonic valence quarks and their super-symmetric partners. 

^ We show the existence of an energy scale below which the fluctuations of the QCD 

^^ , Dirac spectrum are given by a chiral Random Matrix Theory (chRMT) with the 

^\ • global symmetries of the QCD partition function. Physically, for valence quark 

r-f^ ' masses below this scale the partition function is dominated by the zero momen- 

^\1 ' tum modes. In the theory of disordered systems, this energy scale is known as 

f^^ ' the Thouless energy. In the second part of these lectures we discuss chRMT as 

r^v , a schematic model for the QCD partition function at nonzero temperature and 

^y^ , chemical potential. We discuss novel features resulting from the non-Hermiticity 

■"x^ ■ of the Dirac operator. The analysis by Stephanov of the failure of the quenched 

f | ' approximation, the properties of Yang-Lee zeros, as well as the phase diagram of 

Q.' the chRMT partition function are discussed. We argue that a localization transi- 

I tion does not occur in the presence of light quarks. Several results will be derived 

^1h, in full detail. We mention the flavor symmetries of the QCD Dirac operator for 

QJ ■ two colors, the calculation of the valence quark mass dependence of the chiral con- 

i-i^ ' densate and the reduction of the chRMT partition function to the finite volume 

partition function. 



$—( ' 1 Introduction 



The two main phenomena that characterize the low-energy hmit of QCD are 
confinement and chiral symmetry breaking. Because of confinement and the 
spontaneous breaking of chiral symmetry the low-lying states of QCD are given 
by the corresponding Goldstone bosons. At low temperatures we expect that 
the thermodynamics of the QCD partition function iS-Jialj of a gas of (almost) 
massless pions. Mainly because of lattice simulationsEluoQ we know that chiral 
symmetry is restored above a critical temperature of about 140 MeV . The 
order parameter of the chiral phasa transition is the chiral condensate which 
through the Banks-Casher formulau is directly related to the spectral density 
of the QCD Dirac operator. This is our main motivation for detailed studies 
of QCD Dirac spectra. 



One of the questions we wish to address in these lectures is to what extent 
the infrared hmit of the Dirac spectrum can be obtained analytically. After all, 
the infrared limit of the QCD partition function is determined by the chiral 
Lagrangian describing Goldstone modes that become noninteracting at zero 
momenta. On the other hand we know from many studies of complex systems 
that the correlations of eigenvalues on the scale of the individual level spacings 
are universal and are given by Random Matrix Theory. More jprecisely, this 
has been formulated as the Bohigas-Giannoni-Schmit conjectureQ which states 
that spectral correlations of a classically chaotic system are given by Random 
Matrix Theory (see Guhr et al. U for a comprehensive review) . 

The natural question one may ask is whether one can reconcile the uni- 
versal nature of the chiral Lagrangian with the universal spectral correlations 
given by RMT. The first hint for the existence of a close connection came 
from the work of Leutwyler and Smilga a who found that the existence of a 
low-energy chiral Lagrangian imposes an infinite family of constraints on the 
spectrum of the Dirac operator. On the other hand, in the field of Random Ma- 
trix Theory it has been known for a long time that the spectral correlations can 
be described by a Goldstone manifold arising from the spontaneous brrjaking of 
symmetry between the advanced and retarded Greens functions BliEllid'EJ. The 
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through in thia^problem came with the introduction of valence quarks 



|i3||ij.[_5[|i6|£7| r^Yie idea isEa to study a partition function that, in addition to the 
usual sea-quarks, contains a valence quark and its bosonic super-partner. After 
differentiation with respect to the source term and equating the fermionic and 
bosonic valence quark masses, the determinants with the valence quark masses 
cancel and we have access to the Dirac spectrum weighted by the standard 
QCD action. As is the case for the usual chiral Lagrangian, at low energies the 
partition function involving both sea quarks and valence quarks is completely 
determined by the scheme of chiral symmetry breaking and Lorentz-invariance 
and has_heen analyzed in the context of partially quenched chiral perturbation 
theory EjS However, instead of constraints on QCD-Dirac spectra such as 
were obtained in the work by Leutwyler and Smilgau, we are now in a position 
obtain an exact analytical expression for the Dirac spectrum in the domain of 
applicability of this so called partially quenched chiral pertuibation theory. 

As is the case for the usual chiral perturbation theorycil, one can distin- 
guish an important mass-scale. For sufRciently large volumes such that the 
contributions of the non-Goldstone modes are suppressed and valence quark 
masses 



(the pion decay constant is denoted by F, E is the chiral condensate and L 



^^ ^ ^T72 (1) 
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is the length of the box) the contribution of the nonzero momentum modes 
factorizes from the partition function and its mass '^'^Pj!^!]'^!}'"?, Yim^^^*^" ^^ ^^^ 
contribution from the zero momentum modes alone ESEjO'Ejil^'U. However, 
QCD is not the only theory with this zero momentum sector. Any theory with 
the same pattern of spontaneous symmetry breaking as QCD can be reducpij 
this zero momentum limit. In particular, chiral Random Matrix Theory t 
with randomly distributed matrix elements of the Dirac operator is such a 
theory. In this case the strong interactions required to achieve a spontaneous 
breaking of chiral symmetry result from strongly non-local random interac- 
tions. 

Such a picture is well-know from the theorjj^aLiiifiordered mesoscopic sys- 
tems (a number of recent reviews are available LtEjE]). In this case, the scale 
below which Random Matrix Theory is valid is known as the Thoulesa_energy, 
which is the inverse diffusion time of an electron through the sample l3 

E. = ^, (2) 

where D is the diffusion constant and L is the linear dimension of the box. The 
domain with SE <^ Ec is known as the ergodic domain whereas the domain 
Ec ^ SE ^ h/Te (with Te is the elastic scattering time) is called the diffusive 
domain or the Altshuler-Shklovskii domain. In the case of QCD, the diffusion 
process is that of a quark propagating through the Yang-Mills fields in a 4+1 
dimensional space time. The interpretation of spontaneous chiral symmetry 
breaking as a Mott transition was given earlier in£p. By analogy with the 
Kubo formula, E plays the role of the conductivityEj. 



By now, this picture of the QCD gauge field configuratioim 
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g|itimJiipiJ|H.:iprimjcsti£ated in many lattice QCD simulationsliJ'LiufrlEffiila 
@Bm B@0'0@ and several instanton liquid simulations Hb^Hy. In 
particular, it has been confirmed that the scale below wkich_c]|iir|al Random 
Matrix Theory is valid is given by the Thouless energy ([l|)ll30EJ^E3ca. It turns 
out that eigenvalue correlations are given by chRMT all the way up to this 
scale. Beyond this energy scale one has to take into account the contributions 
of the nonzero momentum modes. This can be done simply to one-loop or- 
der in chiral perturbation theory. In that .Wfy, one obtains analytical results 
for the valence quark mass dependence e2H13 of the chiral condensate in the 
Altshuler-Shklovskii domain. 

In the first part of these lectures we will show that chRMT is an exact the- 
ory for the spectrum of the QCD Dirac operator in the ergodic domain. This 
is achieved by showing that in the range (0) the valence quark mass depen- 
dence of the chiral condensate obtained from the low-energy limit of the QCD 



partition function coincides with the valence quark mass dependence obtained 
from chiral Random Matrix Theory. In the second part of these lectures we 
take a diffeijen.t. rQii.te.an.d use chRMT as a schematic model of the chiral phase 
transitionEdOSEJO'Ell. In that case we obtain only qualitative results for the 
spectrum of the QCD Dirac operator and the chiral phase transition. For ex- 
ample, such models are the random matrix equivalent of a Landau-Ginzburg 
functional and give rise to mean field critical exponents. At the tricritical 
point where the upper critical dimension is three they nipvide us with a reli- 
able description of the the infrared degrees of freedom l3. The advantage of 
using random matrix models as opposed to a Landau-Ginzburg functional is 
that one has access to the Dirac spectrum. In particular, at nonzero chem- 
ical potential, this has provided us with valuable insights in thfi-fjSects of 
non-Hermiticity on the properties of the QCD partition function eJ'E2I. The 
common ingredient of chRMT and QCD at nonzero chemical potential is that 
the phase of the fermion determinant due to the non-Hermiticity makes Monte 
Carlo simulations impossible. Since the chRMT partition function can be eval- 
uated analytically, we are in a position to analyze some of these problems in 
great detail. In this way, a complete analytical understanding of the failure of 
Ihe quenched approximation at nonzero chemical potential has been obtained 

In the first lecture (section 2) we summarize the main properties of the 
Euclidean QCD partition function. We discuss in detail the global symrnfcj 
tries and give a brief review of lattice QCD (see lectures by Tom Blum e3 
for more details). In the second lecture we analyze the infrared limit of the 
QCD partition function (section 3). We discuss the domain of validity of the 
zero momentum limit of the effective chiral partition function and introduce 
the partially quenched chiral partition function. An important point of this 
lecture is the discussion of the integration manifold. Explicit calculations of 
the valence quark mass dependence of the chiral condensate that rely on the 
structure of this Riemannian super-manifold are presented in section 4. Re- 
sults are derived in detail both by perturbative and non-perturbative methods. 
In the third lecture we introduce chiral Random Matrix Theory. We discuss 
its main properties and emphasize the importance of universality (see section 
5). In the second half of this lecture we present lattice QCD results for the 
statistical properties of Dirac eigenvalues and compare the results with random 
matrix theory. In the last lecture we introduce chiral Random Matrix Theory 
as a schematic model for the chiral phase transition. In the framework of this 
model, we discuss chiral symmetry breaking, quenching at nonzero chemical 
potential, Yang-Lee zeros in the complex chemical potential plane, the struc- 
ture of Dirac spectra at nonzero chemical potential and the phase diagram 



in the mass, chemical potential and temperature plane. Finally, we present 
general arguments that a localization transition does not occur in QCD with 
light quarks. Concluding remarks are made in section 8. 

2 The QCD Partition Function 

The QCD partition function describing strong interactions in a box of volume 
V3 = L^ can be expressed in terms of the eigenvalues of the QCD Hamiltonian 
E}. as 

ZQCD^Y^^-f3E.^ (3) 



where j3 is the inverse temperature. At low temperatures {j3 -^ 00) the parti- 
tion function is dominated by the lightest states of the theory, namely the vac- 
uum state, with an energy density of -E0/V3, and massless excitations thereof. 
The partition function Zqqd can be rewritten as Euclidean functional integral 
over the nonabelian gauge fields A^ 

ZQCD ^ j dA^,\[ dct{D + mf)e-^^'" , (4) 

where Sym is the Yang-Mills action. We assume that this partition function is 
properly regularized, for example, by a lattice regularization. The gauge fields 
Afj_ are nonabelian gauge fields in Euclidean space time of volume V = L^f3. 
In terms of the generators, T^, of SU{Nc) they can be written as 

a, = a;'^, (5) 

The anti-Hermitean Dirac operator is defined by 

D^j^{d^ + iA^), (6) 

where the 7^ are the Euclidean Dirac matrices with {7^^,71^} ~ 'i.Sf^i,. In these 
lectures we will use the chiral representation of the gamma matrices (with 
diagonal 75). 

2.1 Symmetries of the QCD Partition Function 

It is well-known that the QCD action is greatly constrained by gauge symme- 
try, Euclidean Poincare invariance and renormalizability. These symmetries 



determine the structure of the Dirac operator (0). In this section we will dis- 
cuss the global symmetries of the Dirac operator. They play an essential role 
in its spectral properties in the deepest infrared. In particular, the chiral sym- 
metry, the flavor symmetry and the anti-unitary symmetry of the continuum 
Dirac operator are discussed. 

The chiral symmetry, or the ?7a(1) symmetry, can be expressed in terms 
of the anti-commutation relation 

{l5,D}^0. (7) 

This implies that all nonzero eigenvalues occur in pairs ±iXk with eigenfunc- 
tions given by (f>k and j^^cfik- If Afc = the possibility exists that 75 0^ ~ (t>k, 
so that Afc = is an unpaired eigenvalue. According to the Atiyah-Singer 
theorem, the total number of such zero eigenvalues is a topological invariant, 
i.e., it does not change under continuous transformations of the gauge field 
configuration. Indeed, this possibility is realized by the field of an instanton 
which is a solution of the classical equations of motion. On the other hand, it 
cannot be excluded that A^ = while (j)k and 751/)^ are linearly independent. 
However, this imposes additional constraints on the gauge fields which will be 
violated by infinitesimal deformations. Gcnerically, such situation does not 
occur. 

In a decomposition according to the total number of topological zero 
modes, the QCD partition function can be written as 

zQ^^ = J2e"''Z^^^, (8) 

where 

Z?^^^{l[m'}Y[{Xl+m})).. (9) 

/ k 

Here, (• • ■)i, denotes the average over gauge-field configurations with topolog- 
ical charge i' weighted by the Yang-Mills action. In we introduce complex 
conjugated right-handed and left-handed masses we observe that the 9 depen- 
dence of the QCD partition function is only through the combination me^^^^f . 
A second important global symmetry is the flavor symmetry. This sym- 
metry can be best explained by writing the fermion determinant in the QCD 
partition function as a functional integral over Grassmann variables, 

l[det{D + mf)= /'dV'dVie/''"^"-'^'^''^"^^'^'. (10) 



In a chiral basis with i/j^, = 751/' i?, and V-'l — ~"f5ipL, this can be rewritten as 

(11) 



m/ = 


= we have the symmetry 




^L ^ Ui^L, ^R ^ ^i?t/-\ 




V'L ^ V'lV^~\ V'fl ^ ^V-fl- 



(12) 

The only condition to be imposed on U and V is that their inverse exists. If 
the number of lefthanded modes is equal to the number of right-handed modes 
we thus have an invariance under Gl{Nf)]i x Gl{Nf)i^. However, if the number 
of left-handed modes is not equal to the number of right-handed modes, the 
axial-symmetry group is broken to an Sl{Nf) subgroup whereas the vector 
symmetry with U — V remains unbroken. For m — the flavor symmetry is 
thus broken explicitly to Glv{Nf) x SlA{Nf) by instantons or the anomaly. 

What is much more important, though, is the spontaneous breaking of the 
axial flavor symmetry. From lattice QCD simulations and phenomenological 
arguments we know that the expectation value 

i^^P) = {^Pr4>r) + (VJiVL) ~ (240 MeVf (13) 

in the vacuum state of QCD instead of the symmetric possibility {ipip) = 0. 
Phenomenologically, this is known because of the presence of Goldstone modes. 
The pions are much lighter than the a mesons. The spontaneous breaking of 
the axial symmetry also follows from the absence of parity doublets. For 
example, the pion mass and the 6 (or uq) mass are very different (m^ = 
135MeV and ms = 980MeV)^ 

On easily verifies that (iptp) is only invariant for U — V. The vacuum 
state thus breaks the chiral symmetry down to Glv{Nf). Notice that only 
the axiaLsymmetries are broken. This is in agreement with the Vafa-Witten 
theorem EJ which states that vector symmetries cannot be broken in vector-like 
theories such as QCD. We also observe that tii£ complete axial group is broken. 
The reasons behind this maximum breaking Ej of chiral symmetry are less well 
understood. The Goldstone manifold is given by the maximum Riemannian 
submanifold of SlA{Nf) which is the usual manifold SUA{Nf). The complex 
extension of SU{Nf) does not give rise to additional conserved currents, and 
therefore, the total number of Goldstone modes remains the same (we thank 
J.C. Osborn for this remark). 



For three or more colors the gauge fields are not related by complex con- 
jugation and there are no anti-unitary symmetries. However, the situation is 
more interesting for Nc = 2 which will be discussed in the next section. 

2.2 Special Role of Nc = 2 

The group SU{2) is the only nontrivial special unitary group that is pscudoreal. 
As will be shown in this section, the consequence is that for Nc = 2 the 
symmetry group of the QCD partition function in the chiral limit is enlarged 
to U{2Nf). 

As before, the fermionic action is obtained by writing the determinant in 
the partition function as a Gaussian integral over the Grassmann fields (f> and 
(j). Using the chiral representation of the 7 matrices with 

..^(,^-). (14) 

and a^ = {l,i(7k) with ak the Pauli cr-matrices, the fermionic action can be 
written as 

Nf 



/=i 



L , 

(15) 



For Nc — 2, we have that Aj = — r2A^T2 (with T2 = o"2 in color space). 
Combining this with the relation a* = (J2(^^i.<^2 we find 

^L^id^, + iA^)(t)^j^ = (j){^a2T2af,{df, + iA^,)(j2T2^{, (16) 

where we have used that 9^ is anti-Hermitean and that the fermion fields 
are anti-commuting Grassmann variables. The fermionic action can thus be 
rewritten as 

N f — £ £ 



J t^\cr2T2K 



j^^ ,Cr2T2(pR/ V (7^(9^+iA^)y \(72T2(l)i 



(17) 



Obviously the symmetry group is Gl{2N f). If the number of left-handed modes 
is not equal to the number of right-handed modes because of the anomaly, an 
axial C/(l) is broken explicitly as for three or more colors. 



The mass term is given by 



/=1 



/ \ / n _^,^„-^„\ / ^} 



L \ -m/Cr2T2 



a2T24'iJ \mfa2T2 J \a2T2(t)i 



0^ \ / mfC72T2\ f (t>R 

O-2T20^/ \~mfa2T2 ) Vcr2T20^ 



(18) 



Also in this case we expect m.axinium spontaneous chiral symmetry break- 
ing consistent with the Vafa-Witten theorem. This means that only the sub- 
group of U(2Nf) that leaves both 4'r,4>r and (pLfj^L invariant remains unbroken. 
Therefore, only the subgroup that leaves 



!:6t a2T2 T ) a.r'U ' ^^^^ 



invariant remains unbroken {fpRipR is invariant under the same transforma- 
tions). This is the symplectic group Sp{Nf). The Goldstone manifold is thus 
given by the coset SU {2N f) / Sp{N f) . 

For Nc = 2 the Dirac-operator has the anti-unitary symmetry 

[KCT2,iD]=Q, (20) 

where K is the complex conjugation operator, C — 7274 is the charge conju- 
gation matrix. Because 

{KCT2f = 1 (21) 

we can always find a basis such that the Dirac matrix is real for any D. The 
proof is along the same lines as the proof that time reversal symmetry results 
in real matrix elements for the Hamiltonian in quantum mechanics. 

The anti-unitary symmetry would be violated in the presence of gauge 
fields B^ coupling to the axial current (which are not present in QCD). In that 
case the Dirac operator is given by 

1^L^^, + ilt^Af, + 7^75^^' (22) 

which, because of the absence of the factor i in the _B^-term does not satisfy the 
commutation relation (pOJ). In other words, the invariance of the Dirac operator 
under (kO) follows from the condition ^^D^^ — 2D (in four dimensions). For 
Nc > 3 this condition does not affect the anti-unitary symmetries, and the 
the matrix elements of the Dirac_operator are arbitrary complex numbers both 
with and without the B^, termE3. 



2.3 Dirac Operator in a C'hiral Basis 

It is instructive to write the Dirac operator in a chiral basis with n right- 
handed modes (p^ and n+iy left-handed modes cj)-^. Using the anti-commutation 
relation m) one can easily derive that 

{<pf\D\cj^f) = (0f |75i575|0f ) = -(0f l^l^f ) (23) 

from which it follows that (^^iZ^l^f ) = 0. Similarly, (0^|Z?|0f ) = 0, resulting 
in a Dirac Matrix with the structure 

^-(-Tto)' (24) 

where T is an n x (n + v) matrix. From the inspection of the secular equation 
one concludes that such matrix has always exactly v zero eigenvalues. For 
example, the eigenvalues of the matrix 



(25) 




are equal to and ±iy/a*a + b*b. The matrix structure (e3) is at the basis 
of the Random Matrix Theory to be discussed in section 4. A particular 
realization of a chiral basis is a basis of instanton zero modes whichJias been 
used extensively in the instanton liquid model of the QCD vacuum oE3 

24 Lattice QCD 

In lattice QCD, the gauge fields represented by unitary matrices on the links 
of the lattice, i.e. 

U,^,+, = e^*^^(^), (26) 

where /i denotes a unit vector in the direction of x^ . The gluonic lattice action 
of QCD is the plaquette action given by 

S^IY.(Sp + sU- (27) 

p 

For example, the action of an elementary plaquette 1234 is given by 

Sp = TrC/i2[/23f/43'(7f/. (28) 

10 



This action is invariant under the local gauge transformations 

Uki -> V^^^UkiVi. (29) 

The fermion fields are represented by Grassmann valued fields on the sites. 
Under gauge transformations they transform according to 

V-fc ^ Uki^i. (30) 

However, a lattice discretization of the Dirac action is more problematic. 
In a naive lattice discretization we are faced with the ferpaJon doubling problem. 
This is highlighted by the Nielsen-Ninomiya theoremE3, which states that a 
lattice discretization of chiral fermions with local interactions is not possible. 

There are two widely used discretizations of the fermionic action which 
deal with the doubling problem at the expense of some of the symmetries of 
the continuum QCD action: the Kogut-Susskind or staggered action and the 
Wilson action. The Kogut-Susskind action is defined by 



■^^"^ = n^'n^J.{'n■)Xn{Un,n+^,Xn+^, - Ul_^^^^Xn-f,) + m) ^XnXr, 



n ^ n 

= XmD^nXn, (31) 

where '7^(fi) is a phase factor resulting from the diagonalization of the 7- 
matrices. For m — and ni +n2 + n^ + n^ even, this action is invariant under 
the transformations 

(32) 



Xn -^ e'\n, 


Xn+fj. -^ e Xn+ti 


all rii we have the invariance 




Xn ' c ^„, 


Xn -> e-'"Xn. 



(33) 

For m 7^ this U{1) x U{1) symmetry is broken to the second U{1) symmetry 
only. If we organize the Kogut-Susskind Dirac matrix according to blocks with 
even and odd ni + n2 + n^ + ^4, one obtains the same block structure as in the 
chiral representation of the continuum Dirac operator (p4|). This shows that 
that the eigenvalues of D^^ occur in pairs ±A. 

The staggered Dirac operator is anti-Hermitean with eigenvalues iXk on 
the imaginary axis. The Dirac operator satisfies the sum-rule 

TtD'^'^D^^ = TrD^^D^^I^^^ . (34) 

11 



This follows from the observation that U only occurs in the combination UU^^ 
in the left hand side of this equation. An important application of this sum 
rule is as a test of the accuracy of numerically calculated eigenvalues of the 
Dirac operator. 

The Wilson Dirac action is defined by 

n fj, n 

= ^PiDZ^^i. (35) 

The term proportional to r was introduced by Wilson to remove the fermion 
doublers. However, it also destroys the Ua{^) symmetry of D^ , and the 
eigenvalues do not occur in pairs ±A. On the other hand, this action is invariant 
under the flavor group SUv{Nf). The Wilson Dirac operator satisfies the 
Hermiticity relation 

D^f=j,D'^^, (36) 

from which it follows that 75-0*^ is Hermitean. In the literatui:e.pirac spectra 
of both D^ with eigenvalues scattered in the complex planeE3'L3E3 and 75!?^ 
with eigenvalues_an the real axisEj'LJ have been studied (several reviews ap- 
peared recentlyOO) . Below we will restrict ourselves to the Hermitean Wilson 
Dirac operator. The Wilson Dirac operator satisfies the sum rule 

TtD^^'D^ ^ TtD^''D^\jj^^ . (37) 

Again this can be seen from the observation that U only occurs in the combi- 
nation UU~^ in the l.h.s. of this equation. 

The anti-unitary symmetries of the Wilson Dirac operator are the same as 
of the continuum Dirac operator. For Nc = 2 we thus have 

[CKT2,iD^]=Q. (38) 

For staggered fermions only the eigenvalues of the 7 matrices appear in the 
Dirac operator. Since they are real they do not affect the anti-unitary symme- 
try. For A^c = 2 we thus have 



[Kt2,D^^]=0. (39) 

However, (A't2)^ = —1, and it can be shown that it is always possible to 
find a basis such that the matrix elements of D^ can be organized into real 

12 



quaternions for any gauge field configuration. The eigenvalues of sucli matrix 
are scalar quaternions and are thus doubly degenerate. This can also easily 
be shown from (|39|). Because of this relation, if (j) is an eigenfunction with 
eigenvalue A then Kt24> is also an eigenfunction with eigenvalue A. But 

(0, KT2(b) = {Kt2<J), 4>)* = {{KT2f^, Kt2^) = -(0, Kt2^), (40) 

from which it follows that (f) and Kt2<J) are linearly independent. For a more 
detailed discussion of the symmetries of the staggered Dirac operator we refer 



toumA 

2.5 The Chiral Phase Transition in QCD 

It is expected that chiral symmetry will be restored above some critical tem- 
perature. This is supported both by naive arguments based on counting the 
total number of degrees of freedom and by lattice QCD simulations B. In the 
chiral limit, the order parameter for the chiral phase transition is the chiral 
condensate, {ipip). For QCD with three color and two massless quarks, it van- 
ishes above a critical temperature of about 140 MeV. The nature of the phase 
transition is still under dispute. It is not_yet clear whether the transition is 
a second order or a weak first order^oiifiEj. The role of instanton field con- 
figurations has to be clarified as wellLjO. However, with current progress in 
computational resources, these questions should be answered in the near fu- 
ture. The only remaining fundamental problem is the study of QCD at finite 
baryon density. Because of the phase of the fermion determinant, this problem 
cannot be resolved by means of Monte-Carlo simulations. Progress in this area 
requires new paradigms which make it a particularly challenging area of active 
research. 

2.6 The Banks- Casher Relation and Microscopic Spectral Density 

The order parameter of the chiral phase transition, (ipip), is nonzero only below 
the critical temperature or a critical chemical potential. As was shown by 
Banks and CasherH, {"ipip) is directly related to the eigenvalue density of the 
QCD Dirac operator per unit four-volume 

S-IW)l=lim^, (41) 

where the spectral density of the Dirac operator with eigenvalues {Afc} is de- 
fined by 

p(A)^(^5(A-Afc)). (42) 

k 
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It is elementary to derive this relation. The chiral condensate is defined as 
the logarithmic derivative of the the partition function (0). For equal quark 
masses we obtain, 

W> = -lim-i-9™logZQ^^(™) ^ -liml(^-^^). (43) 
VNf V '^ Xj. + m^ 

If we express the sum as an integral over the average spectral density, and take 
the thermodynamic limit before the chiral limit, so that many eigenvalues are 
less than m, we recover mW). The order of the limits in mw is important. 
First we take the thermodynamic limit, next the chiral limit and, finally, the 
field theory limit. As follows from (p3), the sign of {"tp^) changes if m crosses 
the real axis. 

An important consequence of the Bank-Casher formula ( pl| ) is that the 
eigenvalues near zero virtuality are spaced as 

AA = l/p(0) = tt/SF. (44) 

This should be contrasted with the eigenvalue spectrum of the non-interacting 
Dirac operator. Then the eigenvalues are those of a free Dirac particle in a 
box with eigenvalue spacing equal to AA ^ XjV'^l^ for the eigenvalues near 
A = 0. Clearly, the presence of gauge fields leads to a strong modification of 
the spectrum near zero virtuality. Strong interactions result in the coupling of 
many degrees of freedom leading to extended states and correlated eigenvalues. 
Because the eigenvalues near zero are spaced as ~ XJYjV it is natural to 
introduce the microscopic variable 

u = AV^S, (45) 

and the microscopic spectral densitycil 

P.(u) = ^lnn^^(p(^)). (46) 

We expect that this limit exists and converges to a universal function which 
is determined by the global symmetries of the QCD Dirac operator. The cal- 
culation of this miivcrsal function from QCD is the main objective of these 
lectures. We will calculate Ps(u) both from the simplest theory in this uni- 
versality class which is chiral Random Matrix Theory (chRMT) and from the 
partial quenched chiral Lagrangian which describes the low-energy limit of 
the QCD partition function. We find that the two results coincide below the 
Thouless energy. 
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2. 7 Valence Quark Mass Dependence of the Chiral Condensate 

Instead of studying the Dirac spectrum it is often convenient to consider the 
valence quark mass dependence of the chiral caadepaate. In terms of the 
eigenvalues of the Dirac operator it is defined by E3E3IH3 



S(TO„;mi,---,mjv,) = ^ I] ( 77— 



V ^ \ iXk + ruy 

k 

^^C^ p{X-.m,---,m,,) ^ (47) 

Vj i\ + m^ ' ' 

Here, (• • •) denotes an average with respect to the distribution of the eigenval- 
ues. 

The relation ( ^ ) can then be inverted to give p{X;mi,. . . ,mNj). As 
mentioned inll3, the spectral density follows from the discontinuity across the 
imaginary axis, 

Discl^^^,^ I](m„) = lim S(zA + e) - S(iA - e) = 2t:^{5{\ + A^)) = 27rp(A), 

^^ k 

(48) 
where we have suppressed the dependence on the sea-quark masses. 

3 Infrared Limit of the QCD Partition Function 

3.1 The Chiral Lagrangian 

For light quarks the low energy limit of QCD is well understood. It is given by 
the chiral Lagrangian describing the interactions of the pscudoscalar mesons. 
The reason is that pions are Goldstone bosons which are the only light degrees 
of freedom in a confining theory such as QCD. To lowest order in the quark 
masses and the momenta, the chiral Lagrangian is completely dictated by chiral 
symmetry and Lorentz invariance. In the case of Nf light quarks with chiral 
symmetry breaking according to SUl{Mj) x SUR{Nf) -^ SUv{Nf) the so 
called Weinberg Lagrangian is given byO 

Coff{U) = ^ Tvid^U 9^C/t) - I Tt{MU + MU^), (49) 

where F is the pion decay constant, E is the chiral condensate and M is the 
quark mass matrix. The fields U{x) are SU{Nf) matrices parametrized as 

U = exp(iV2nai7F), (50) 

15 



with the generators of SU{Nf) normalized according to Trt°'t'' — (5°''. This 
chiral Lagrangian has been used extensively for the description of pion-pion 
scattering amplitudes. 

To lowest order in the pion fields the chiral Lagrangian can be expanded 
as (for equal quark masses) 



^d.u^d^n^ + ^- 



J^csiu) = -a^n'^a^n'^ + ^^^'^^^ (51) 



This results in the pion propagator l/(p^ + m^) with pion masses given by the 
Gellmann-Oakes-Renner relation 

^l^^^- (52) 

In the long-wavelength limit the order of magnitude of the-different terms 
contributing to the action of the chiral Lagrangian is given byl^ 

S = l d'xC{U) ^ L'^-'UlzM^ + i'f?(n|M' + n%zM\ (53) 

Here, the H^j^j represent the zero momentum modes with no space time depen- 
dence, whereas the nonzero momentum modes are denoted by by n^2j(^(a;). 
This decomposition has two immediate consequences. First, for ^?^ ^ ^ the 
fluctuations of the pion fields are small and it is justified to expand U in powers 
of n'^. Second, for 

Sto 1 

the fluctuations of the zero modes dominate the fluctuations of the nonzero 
modes, and only the contribution from the zero modes has to be taken into 
account for the calculation of an observable. In this limit the so called finite 
volume partition function is given by 

ZfAM,0)^ f dC/e^^^^T'^•^^^"''"^ (55) 

Juesu(Nf) 

where the ^-dependence follow^ from the (?— dependence of the QCD partition 
function via the combination^ me^^^^f . We emphasize that any theory with 
the same pattern of chiral symmetry breaking as QCD can be reduced to the 
same extreme infrared limit. 
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3.2 Leutwyler-Smilga Sum Rules 

The finite volume partition function in the sector of topological charge v follows 
by Fourier inversion according to (^ . The partition function for ly — is thus 
given by ( p5[) with the integration over SU{Nf) replaced by an integral over 

The Leutwyler-Smilga sum-rules B are obtained by expanding the parti- 
tion function in powers of m before and after averaging over the gauge field 
configurations and equating the coefficients. This corresponds to an expansion 
in powers of tti of both the QCD partition function (Q) and the finite volume 
partition function ( [55| ) in the sector of topological charge z/. As an example, 
we consider the coefficients of m^ in the sector with i/ = 0. After performing 
the group integrals we find the sum-rule 



(E 






(56) 



where the prime indicates that the sum is restricted to nonzero positive eigen- 
values. 

By equating higher powers of m^ one can generate an infinite family of 
sum-rules for the eigenvalues of the Dirac operator. However, they are not 
sufficient to determine the Dirac spectrum. The reason is that the mass in 
the propagator also occurs in the fermion-determinant of the QCD partition 
function. However, as will be shown in the next section, the Dirac spectrum can 
be obtained from a chiral Laepangian corresponding to QCD with additional 
bosonic and fermionic quarks Ii3. We conclude that chiral symmetry breaking 
leads to correlations of the inverse eigenvalues which are determined by the 
underlying global symmetries. 

3.3 The Partially Quenched QCD Partition Function 

As was mentioned in previous section, the valence quark mass dependence of 
the chiral condensate cannot be extracted from the QCD partition function. 
The solution to this problem i^o simultaneously introduce yet another quark 
species of opposite statistics t£l. This corresponds to a Euclidean partition 
function of the form 



Nf 

^p, = fdA t"lin^""''\ n det(i? + -/) e"^-^ , (57) 

J det{D + m^2) f}^ ^ 
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which we wih cah the partiaUy quenched or pq-QCD partition function. When 
rUyi = mv2 this partition function simply coincides with the original QCD 
partition function. However, it is now also the generator of a mass-dependent 
chiral condensate (see (^7|)) for the additional (say, fermionic) quark species. 
In the sector of topological charge v we find 



E(m^;mi,---,mAr^) 



d 



1 



V dniyi 



log^r, 



(58) 



rrii,!— 771^,2— rrii, 



where ZP^^ is the partially quenched QCD partition function in the sector of 
topological charge v. 

Our aim is to find the chiral Lagrangian corresponding to (p7h. If we are 
successful, we have succeeded in deriving a partition function for the extreme 
infrared limit of the spectrum of the QCD Dirac operator. These questions 
will be addressed in the next sections. 



3.4 The Infrared Limit of QCD 

If we don't write the determinants in terms of integrals over complex conju- 
gated variables, the global flavor symmetry of the partially quenched QCD 
partition function (|57|) is broken spontaneously according to 



(59) 



GlniNf + 1|1) X GhiNf + 1|1) ^ Glv{Nf + 1|1) 



with an axial U{1) group broken explicitly by the anomaly. Here, the groups 
GlniNf + 1|1), GhiNf + 1|1) and Glv{Nf -I- 1|1) are super-groups of matrices 
acting on vectors with Nf + 1 fermionic components and one bosonic com- 
ponent. The subscript refers to right-handed (R), left-handed (L) and vector 
(V), respectively. The latter transformations transform the right-handed and 
the left-handed fermion fields in the same way. For a confining theory such as 
QCD the only low-lying modes are the Goldstone modes associated with the 
spontaneous breaking of chiral symmetry. The quark masses play the role of 
symmetry breaking fields. 

Although the axial supergroup GlA{Nf -I- 1|1) is a symmetry group of the 
pq-QCD action ( pT] ) it is not necessarily a symmetry of the QCD partition 
function. It may be that the symmetry transformations violate the conver- 
gence of the integrals in the partition function. There are no problems for the 
Grassmann variables. However, the integrations over the bosonic fields are only 
convergent if the fields are related by complex conjugation. For this reason a 
[/a(1) transformation in which the bosonic fields 4>r and (pL are transformed 
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according to a different phase factor is not a symmetry of the partition func- 
tion. What is a symmetry of the partition fmiction is the axial transformation 

0L -* e^cjiL, (j)*L -^ e''(j)l, (60) 

^R ^ e-'^R (^R ^ e-'cj)R. (61) 

Mathematically, this symmetry group is Gl{\)/U{\). Had we restricted our- 
selves to the unitary subgroup U{Nf -I- 1|1) of Gl{Nf + 1|1) from the start, we 
would have missed this class of symmetry transformations. 

Taking into account the chiral anomaly, the chiral symmetry in (m) is 
broken spontaneously according to (^9|). The symmetry of the QCD parti- 
tion function is thus reduced to Slv{Nf -I- 1|1) Glv{i) where denotes the 
semi-direct product. The Goldstone manifold corresponding to the symmetry 
breaking pattern ( p9|) is based on the symmetric superspace SlAiNf + 1|1). 
In our effective partition function the terms that break the axial symme- 
try will be included explicitly resulting in an integration manifold given by 
SlA{Nf + 1|1) GIa{^)- However, this manifold is not a super-Riemannian 
manifold and is not suitable as an integration domain for the low energy parti- 
tion function. As an integration domain we choose the maximum Riemannian 
submanifold of GlA{Nf + 1|1). This results in a fermion-fermion block given 
by the compact domain UA{Nf + 1), whereas the boson-boson block is re- 
stricted to the non-compact domain GI{1)/U{1). Because of the super-trace, 
this compact/non-compact structure is required for obtaining a positive defi- 
nite quadratic form for the Jiiass term and the kinetic term of our low energy 
effective partition functioned. For a detailed piathematical discussion of this 
construction we refer to a paper by ZirnbauerS. 

To lowest order in the momenta, the infrared limit of the QCD partition 
function is uniquely determined by the geometry of the Goldstone manifold 
and Lorentz invariance. With the singlet field $o = Str<I> fluctuating about 
the vacuum angle with an ampl^fude given by the singlet mass we thus obtain 
the effective partition function 



,M) = I dC/exp / d'^x 

'ueGl{Nf + l\l) 



^ str{d^u d^u-^) 



+ fstr(^^ + ^t/-) + ^(^ 



(62) 



This partition function contains terms up to order p^ . The kinetic term of the 
singlet field is subleading and has not been included. The supertrace (Str) and 

19 



the superdeterminant (Sdct) of a graded rnatdx with bosonic blocks a and b 
and fermionic blocks a and p arc defined as 



IS 



Str ( "" ^ ) = Tra - Tr6, (63) 



and 



det(a — ab ^ p) 

Sdet = ^ — ; ^. (64 

det 6 ^ ' 

Recently, this partition function for N f — was derived from a two sub- 
lattice randoraJlux model E3 using the flavor-color transformation introduced 
by Zirnbauer eZI. In our language, such model is QCD at infinite coupling 
with U{N) gauge fields and quenched Kogut-Susskind fermions with no phase 
factor from the 7 matrices. In this case, as well as in the effective partition 
function for Kogut-Susskind fcrmion, the singlet mass term is absent. In both 
cases the f/(l) x f/(l) symmetry is broken spontaneously to U{1) instead of 
broken explicitly by the anomaly. The reason is that the U{1) symmetry of 
Kogut-Susskind fermions is an axial isospin transformation rather than the 
singlet Ua{^) transformation (We thank Paul Rakow and Misha Stephanov for 
discussions to clarify this point.). I— . 

By a generalization of an argument due to Gasser and LeutwylerEil to be 
discussed in the next section it can be shown that for masses in the range (El) 
this partition function factorizes iijitp a product over zero momentum modes 
and non-zero momentum modes tSll3. To leading order, the mass dependence 
of the of the QCD partition function is thus given by 

JueGl{Nf + l\l) 

The partition function in a sector of topological charge v follows by Fourier 
inversionB 

1 c27r 



ZesiM)^— dee-''''Z,s{e,M). (66) 

We perform the integration over 6 after a Hubbard-Stratonovitch transforma- 
tion which linearizes the singlet masajterm. Up to a mass independent factor 
this results in the partition function t3 

Z^eiM)^ f dUSdefUe^'istriMU+MU-')^ (g7) 

JueGl{Nf + l\l) 
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3.5 Domains in (Partially- Quenched) C'hiral Perturbation Theory 

In chiral perturbation theory the different domains of vahdity where analyzed 
in detail by Gasser and Leutwyler Eil. A similar analysis applies to partially 
quenched chiral peLturbation theory. The idea is as follows. The U field can 
be decomposed asEJ 

U = Uoe'^'^'=\ (68) 

where Uq is a constant (zero-momentum) field. The kinetic term of the low- 
momentum components is the f/' fields can be approximated by 

]^jd^xd,,rd^r-L^^^- (69) 

We observe that the magnitude of the fluctuations of the ip field are of order 
1/ L which justifies a perturbative expansion of exp(j'0(a;)). The fluctuations 
of the zero modes (i.e. constant fields), on the other hand, are only limited by 
the mass term 

]-VT.StvM{Uo + U^^). (70) 

For quark masses m ^ l/VS, the field t/o fluctuates close to the identity and 
the C/p field can be expanded around the identity as well. This is the domain 
of chiral perturbation theory. 

For valence quark masses in the range 

^ « m. « ^ (71) 

the valence quark mass dependence of the chiral condensate is dominated by 
the zero momentum modes which can be treated perturbatively. Below we will 
show that in this domain chiral perturbation theory and random matrix theory 
coincide. In the theory of disordered mesoscopic systems it is well-known that 
random matrix theory is valid below an energy scale ~ 1/L^. This will be 
discussed in the next section. 



3.6 Picture from Mesoscopic Physics 

In disordered mesoscopic physics it has been found that Random Matrix The- 
ory is valid for eigenvalues separated less than the energy scale Ec (see eq. 
(^). In this context Ec is defined as the inverse tunneling time of an electron 
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through the sample. For diffusive motion, the distance Ax a particle diffuses 
in time interval At is given by 

(Ax^ = DAt, (72) 

where D is the diffusion constant. The diffusion time through the sample is 
thus given by L^/D resulting in a Thouless energy of 

E. - §. (73) 

A second energy scale that enters in mesoscopic physics is h/r,,, where 
Te is the elastic scattering time. Based one these two scales three different 
domains for the energy difference, SE, that enters in the two-point correlation 
function, can be distinguished: the ergodic domain, the diffusive domain or 
Altshuler-Shklovskii domain and the ballistic domain. For energy differences 
6E <C Ec eigenvalue correlations are given by RMT. Since for time scales 
larger than the diffusion time, an initially localized wave packet explores the 
complete phase space, this regime is known as the ergodic regime. In the 
diffusive regime defined by Ec ^ SE ^ h/re only part of the phase space is 
explored by an initially localized wave packet, resulting in the disappearance 
of eigenvalue correlations. In this paper we don't consider the ballistic regime 
with 6E ^ Ti/tc- For an interpretation of the Thnuleas energy in terms of the 
spreading width we refer toBEj. As was shown inE3'Eil, the spectral two-point 
function can be related to the semiclassical return probability which provides a 
simple intuitive pirtiirp nf its a^V'^r'^p^i'" behavior. For other recent studies on 
this topic we refer toL3o'oE3€3'EZlE3. What has emerged from these studies is 
that there is a close relation between eigenvalue correlations and localization 
properties of the wave functions. 

Based on these ideas we can interpret the Dirac spectrum as the energy 
levels of a system in 4 Euclidean dimensions and one artificial time dimension. 
The corresponding classical evolution can be easily derived from the Heisenberg 
equations of motion 

^=z[x,i?]. (74) 

According to the Bohigas conjecture the eigenvalue correlations are given by 
RMT if and only if the corresponding classical motion chaotic. We thus con- 
clude that the classical time evolution of quarks in the Yang-Mills gauge fields 
is chaotic. 

These ideas can be tested by means p£. lattice QCD EJC3 and instanton 
liquid Ej simulations. In Fig. 1, we show til the number variance S2(n), de- 
fined as the variance of the number of levels in an interval containing n level 
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Figure 1: The number variance S2(n) versus n approximation for an interval starting at 
A = 0. The total number of instantons is denoted by N. 



on average, versus n for eigenvalues obtained from the Dirac operator in the 
background of instanton liquid gauge field configurations. The chRMT result, 
given by the solid curve, is reproduced up to about two level spacings. In 
units of the average level spacing, A = l//5(0) = irfSV, the energy E^ is given 
by ric = Ec/A — F'^L'^/tt. For an instanton liquid with instanton density 
N/V = 1 we find that n^ « O.OTVN. We conclude that chRMT appears to 
describe the eigenvalue correlations up to the predicted scale. 

In mesoscopic physics the Goldstone modes that enter in the theory of im- 
purity scattering are known as the diffusons. The reason is that the Goldstone 
Sopagator satisfies a diffusion equation (see for example the bonk by Efetov 
). It describes the diffusion of electrons by impurity scattering 113. Similarly, 
the pion propagator also satisfies a diffusion equation. This can be seen by 
Fourier transforming the propagator with regards to the imaginary valence 
masses my = iX and m^/ = i\', 



n(x,T) 



d^p 



dXdX'- 



^ipx~\-iT{X--\-X') 



{2t:Y J p^ + {iX + iX')^/F^' 

For (x, r) 7^ this propagator satisfies the diffusion equation 

{Ddl-dr)n{x,T) =0 
with a diffusion constant given bytila D — F^fS . 
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(75) 



(76) 



4 Calculation of Valence Quark Mass Dependence of the Chiral 
Condensate from QCD 

4-.1 Nonperturbative Evaluation ofE^rriy) 

In this section we calculate the valence quark mass dependence of the chiral 
condensate for the simplest case of Nf = and i/ = in the domain m^, <C 
F"^ /TiL/^ . In this domain the partition function is given by 

Z{J)^ I df/exp[^Str(''"^ + '^ ]{U + U~') 



ueGi(i\i) 







(77) 



where the integration is over the maximum super-Riemannian sub-manifold of 
GZ(1|1). This manifold is parametrized by 

' e"*' 



U = exp 



a 

(3 



e' 



(78) 



The integration measure is the Haar measure which in terms of this parametriza- 
tion it is given by 

5Uki 



Sdet- 



dadl3d(j)ds, 



(79) 



5(j) 6s da 5(3 

where 5U = U-^dU. 

It is straightforward to calculate the Berezinian going from the variables 
dUii SU22 SU12 6U21 to the variables dcj) 6s da 5j3. The derivative matrix is given 
by 

/ i (3/2 a/2 \ 



B = 



SUh 



64>6s5a6l3 



(80) 
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1 (3/2 a/2 
e"-"^ 
\0 6-'*+ 

Using the definition of the graded determinant one simply finds that SdetS = i. 
Up to a constant, the integration measure is thus givenjDy d<j)dsdad(3. In 
general, for Nf ^ 0, the Berezinian is more complicated oO. 
We also need 

(1 + ^) 

/3(e*'^-e-^) coshs(l- 

After differentiating with respect to the source term {Yj{my) — dj log Z{J)\j=o/V) 
this results in 

Y.{my) f d(f>dsdad() 



liu + u-') 



2 J 



(81) 



2tt 



. C0S(?i(l + ^)g^cos0(l+^)-xcoshs(l-^)^ 



(82) 



24 



With the Grassmann integral given by the coefficient of a/3 we obtain 
S(r7i„) f dsdcj) , 



S J 4. ^^^^"^ 

+ a;(cos0 + cosh s) cos 0]e^('=°'"^-'=°''^'). (83) 

Now all integrals can be expressed in terms of modified Bessel functions. We 
find 

^^ = h{x)Ko{x) + ^{l2{x)Ko{x) + Ioix)Ko{x) + 2h{x)K^{x)). (84) 

After using the recursion relation for modified Bessel functions, 

l2{x)^Io{x)~-h{x), (85) 

X 



we arrive at the final resultlHS'tSIl 

^^ = x{I^{x)K^{x) + h{x)K^{x)), (86) 

where x = mVT,. 

This calculation can be generalized to arbitrary Nf and arbitrary v. The 
calculation for arbitrary Nf is much more complicated, but with a natural gen- 
eralization of the factoijized paranietrization, and using some known integrals 



over the unitary groupt£3, one arrives at an expression in terms of modified 
Bessel functions 

^ - x{Ia{x)Ka{x) + Ia+l{x)Ka-l{x)), (87) 

where a — Nf + \i^\- This result is in complete agreement EH with chRMT to 
be discussed below. 

4-2 Perturbative Calculation 

The valence quark mass dependence of the chiral condensate follows from (see 
eq. (H)) 

^im,)^^dj\ogZ^siM)\^^^. (88) 

In the quenched limit this results in (see eq. ( p2|) for the definition S{U{x))) 

S(m„)=/" dU{x)^djStrMiUix) + U-\x))e'^^"'-''^\ (89) 
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Perturbatively to one loop order we can expand U around the identity, and do 
the usual Gaussian integrations for a fiat measure. The trace of the propagator 
of the Goldstone bosons containing valence quarks is given by 

p 

where L is 1 for a bosonic pion and 1 for a ferniionic pion. The valence pion 
mass, Mjr, is given by the Gell-Mann-Oakes-Renner relation 

2m„i; (m„+ms)i; 2m^E 

^™ = — p^> or M^s = -^ , or M,s = — p^"- (91) 

The calculation of the valence quark-mass dependence of the chiral con- 
densate using the partially quenched chiral Lagrangian is basically straightfor- 
ward. The only complication is the "77"' -mass term (the terms proportional to 
TOq in the Lagrangian). However, the quadratic form in the pion fields can be 
diagonalizedllj resulting in an analytical expression for the propagator. In the 
limit that_Lhc 77' mass is much larger than the mass of the Goldstone bosons, 
we findEEiy 

^^ = 1 - ^U]A{Ml) - A{Ml) + [Ml Ml) aM^„ A(Mi)|. 

(92) 
where the trace of the finite volume pion propagator is given by 

p 

Here, the sum is over momenta in a box of volume i^. From a careful analysis 
of the diverging momentum summations it follows that this propagator can be 
expanded in powers of M^ 

where /3i — 0.140,461 for a 4-dimensional hyper-cubic box, and A is the mo- 



mentum cutoff LLHil. In agreement with a naive inspection of the momentum 



sum, the propagator is dominated by the zero momentum term if 

M«i, (95) 
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or the Compton wavelength of the valence pion is much larger than the size of 
the box. As discussed in section (3.1) we thus find that for 

m„ « ^^ (96) 

the valence quark mass dependence of the chiral condensate is described by the 
zero momentum contribution to the trace of the propagator. Let us investigate 
the zero momentum limit of ( |92[) given by replacing the propagator by the first 
term in (|9J). In the chiral limit with rus — > we then find 



E(l 



Nf 



VMl,F^ ■ 



The nonperturbative result for Y,{m^) obtained for a fixed value of the topo- 
logical charge (p7|) v instead of at 6* = can also be expanded in powers of 

^-l-^^^^^ + M + 0(i,), (98) 
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where the term \v\/x results from the factor nv'^^ in the partition function 
We find that to leading order in 1/x the valence quark mass dependence of the 
chiral condensate does not depend on the topological charge and agrees with 
chiral perturbation theory. 

The result for the valence quark mass dependence of the chiral condensate 
(P2) is valid in the diffusive domain as well. In thai_jcase it is possible to 



rederive and to extend results by Smilga and Stern 1123 for the slope of the 
Dirac spectrum. 

4.3 Two-Point Correlation Function 

Let us do a perturbative calculation of the two-point function in the domain 
that is dominated by the zero-momentum modes, i.e. for valence quark masses 
well below F'^ jYil? . The generating function is defined by 

Z{J,J')= f dUe^^^\ (99) 

JU<^Gl(2\2) 
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where the action is given by 



S{U) 



EV 



/m + J 







Str 



\ 



m' + J' 
m 

\ m' J 



(u + u-'). 



(100) 



The^pectral two-point correlation function follows from the partition function 
byN 



{p{X)p{X')) ^ ^Biscdjdj, ^(J, J')l,7=o,j'=o.™=a,™'=.A' 



(101) 



We are interested in the discontinuity of the connected part of the cor- 
relation function. The only contribution is from mesons containing a va- 
lence quark with mass m and a valence quark with mass m' . If we expand 
U = exp i\/2IlaT'^ / F in powers of IlaT°'/F wc find that the connected spec- 
tral two-point function 



5]2 I 

(p(A)p(A'))c = j^^ Disc|„^^^_„,^,^, j-^ 



(102) 



has been expressed in terms of the discontinuity of the valence pion suscepti- 
bility. To obtain this result we have used among others that 



Str 



/1000\ 
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(103) 



and the super-symmetry of the partition function. To calculate the disconti- 
nuity we wish to remind the reader that 



M„ 



{\m\ + \m'\)j: 



F2 



where I^tiI = y^Jm?). Wc finally find 
(p(A)p(A'))c = ^ 



27r2 V(A + A')2 (A-A')2 



(104) 



(105) 



|W±iich agrees with the asymptotic result for the spectral correlation function 
llHJ for the chiral Random Matrix Theory to be discussed next. 
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5 Chiral Random Matrix Theory 

5. 1 Definition of the Model 

The chiral random matrix partition function withJiie global symmetries of the 
QCD partition function as input is defined byOcS 



ZpiM) = / DW Y[ det(I? + nif) 
f=i 



where 



e-^^^'^w'^^ (106) 



^=(^M^to ' (107) 



iW 



and W is a n X m matrix with v = \n — m\ and N — n + m. As is the case in 
QCD, we assume that the equivalent of the topological charge i^ does not ex- 
ceed VjV, so that, to a good approximation, n = N/2. Then the parameter S 
can be identified as the chiral condensate and N as the dimensionless volume of 
space time (Our units are defined such that the density of the modes N/V — 1). 
The matrix elements of W are either real (/3 = 1, chiral Gaussian Orthogo- 
nal Ensemble (chGOE)), complex (/3 = 2, chiral Gaussian Unitary Ensemble 
(chGUE)), or quaternion real (/3 — 4, chiral Gaussian Symplectic Ensemble 
(chGSE)). For QCD with three or more colors and quarks in the fundamental 
representation the matrix elements of the Dirac operator are complex and we 
have /3 = 2. For Nc = 2 and quarks in the fundamental representation the sit- 
uation is more interesting. As discussed in section (2.4) it is always possible to 
find a basis for which the matrix elements of the continuum Dirac operator are 
real for all gauge fields. Then the Dyson index of the corresponding chRMT is 
/3 = 1. For Kogut-Susskind fermions in this case the anti- unitary symmetry is 
different (see section (2.4)) and the matrix elements of the Dirac operator can 
be organized into real quaternions corresponding to a chRMT with /3 — 4. For 
gauge fields in the adjoint representation the gauge fields are real resulting in 
an anti-unitary symmetry also corresponding to the class /? = 4 for any value 

of A^c. 

The reason for choosing a Gaussian distribution of the matrix elements 
is its mathematically simplicity. This model can be generalized to arbitrary 
potential 

TtW^W -> TtV{W^W) (108) 

where V{x) is such that the resulting probability distribution is well defined. 
In section (5.4) we will argue that the interesting properties of chRMT do not 
depend on the choice of V{x). 
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Below we will also discuss the invariant or Dyson- Wigner Random Matrix 
Ensembles. They are defined as ensembles of Hermitean matrices {H} with 
independently Gaussian distributed matrix elements, i.e. with probability dis- 
tribution given by 

Depending on the anti-unitary symmetry, the matrix elements are real, com- 
plex or quaternion real. They are called the Gaussian Orthogonal Ensemble 
(GOE), the Gaussian Unitary Ensemble (GUE) and the Gaussian Symplectic 
Ensemble (GSE), respectively. Each ensemble is characterized by its Dyson 
index P which is defined as the number of independent variables per matrix 
element. For the GOE, GUE and the GSE we thus have f3 = 1, 2 and 4, 
respectively. 

Together with the Wigner-Dyson ensembles and the superconducting ran- 



dom matrix ensembles E££l the chiral ensembles canJpe classified according to 



the Cartan classification of large symmetric spacesEj. 

5.2 Symmetries 

From the structure of the determinant, it follows immediately that the chi- 
ral random matrix partition function has the same global flavor and Ua{^) 
symmetries as the QCD partition function. 

It can be shown that in the domain iwt) the random matrix partition-func- 
tion can be mapped onto the effective finite volume partition function c3. We 
will discuss such derivation in section (7.3) where we extend this model to 
finite temperature. The most detailed studies to this model jaaj;.rp..porfprmed in 
the quenched limit by means of the supersymmetric mcthodL30L3G3. In that 
case the chRMT partition function can be mapped onto a super-symmetric 
nonlinear cr-model. 

In this model chiral symmetry is broken spontaneously with chiral con- 
densate given by E — limAr^oo 7i'/o(A — > 0)/N, where N is interpreted as the 
(dimensionless) volume of space time. For complex matrix elements {(3 — 2), 
which is appropriate for QCD with three Qnjnore colors and fundamental 
fermions, the symmetry breaking pattern isM SU{Nf) x SU{Nf)/SU{Nf). 
For (3 = 1 and 4 the symni^try breaking pattern i&^U{2Nf)/ Sp{Nf) and 
SU{Nf)/0{Nf) respectively^, the same as in QCDE3. 

Finally, one of the reasons for the mathematical simplicity of this model 
is that it has more symmetry than the QCD partition function The chRMT 
partition function is invariant with respect to the unitary transformation 

W^UWV-'^. (110) 
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Since an arbitrary complex matrix can be decomposed as 

W = UKV-^, (111) 

with A a diagonal matrix and U and V unitary matrices, the chRMT par- 
tition function for arbitrary potential ^(A) can be expressed in terms of the 
eigenvalues Afc as 

zf,{M) = I dWM^xDf n Ar^-^ n -/ w^i + -?)e-^ ^^ ^^'^^\ 

k f f,k 

(112) 
where the Vandermonde determinant is defined by 

Ai\l) = l[iXl-\f). (113) 

k<l 

This result greatly simplifies the mathematical treatment of the chRMT par- 
tition function. For example, it allows us to use the orthogonal polynomial 
method. We wish to point out that for /3 = 2 and nif = the ratio of the 
partition function and m^f''^\ depends only on v through the combination 
Nf + ji'l. This duality between flavau-aiid. topology can be understood more 
directly via Itzykson-Zuber integralsll2Zl't22l. 

5.3 Properties of the Random Matrix Model 



The average spectral density that can be derived from (106) has the familiar 



semi-circular shape. As can be easily derived by means of the orthogonal 



po 

by 



, method, the microscopic spectral density for the chGUE is given 

109|,|104|,|25| 

Psiu) - ^ (JHu) - Ja+l{u)Ja^l{u)) , (114) 



where a — Nf + |i^|. The valence quark mass dependence of the chiral conden- 
sate follows from integration over the microscopic spectral density. If we use 



the microscopic variable u = XVTi as new integration variable, equation (47) 
can be rewritten as 

S(mi,) f°° 2udu 1 u 

In the thermodynamic limit, the spectral density can be replaced by the micro- 



scopic spectral density ( |114| ). The integrals over the Bessel functions are known 
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and result in the foUowu 
of the chiral condensatet 



expression for the valence quark mass dependence 



S(a;) 



= x{Ia{x)Ka{x) + Ia+l{x)Ka-l{x)), 



(116) 



where a = Nf + \iy\, x — myVS and /^ and Ka are modified Bessel functions. 
This result is in perfect agreement with the result obtained from the partially 
quenched chiral Lagrangian given in eq. (^7|). The spectraJLcorrelations near 
A = can also be expressed in terms of Bessel functions t2j, whereas in the 



bulk of tViF' spectrum they are given by the invariant random matrix ensem- 
bles 



no 
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The microscopic spectral density and the microscopic spectral 
correlations can also be derived for /? = 1 and /3 = 4 but these two cases are 
mathematically much more complicated. Many other properties of the chRMT 
partition function have been calculated.. Among others, we mention the distri- 



bution of the smallest eigenvalue 11131113 
in the presence of nonzero quark masses 
the chRMT partition function we refer to 



Lie, 
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■oscopic spectral density 
For more discussion of 



5.4 Universality 

The aim of universality studies is to identify observables that are stable against 
deformations of the random matrix ensemble. Not all observables have the 
same degree of universality. For example, a semicircular average spectral den- 
sity is found for random matrix ensembles with independently distributed ma- 
trix elements with a finite variance. However, this spectral shape does not 
occur in nature, and it is thus not surprising that it is only found in a rather 
narrow class of random matrix ensembles. What is surprising is that the 



microscopic spectral density and the microscopic sp 
hlfi Jiiith respect to a much larger class of deformations 



125 



rtra1 cnrrpl 



117 



118 



11£ 



120 



nr^ arp h; 



121 



122 



126 



^ «ta- 
124124 



Two different types of deformations have been considered, those that 
maintain the unitary invariance of the partition functions and those that break 
the unitary invariance. 

In the first class, the Gaussian probability distribution is replaced by 
P{W) ~ exp(— iV^^^ afeTr(W^'l'li!^'^). For a potential with only ai and 02 
different from zero it was shown r^^ that the microscopic spectral density is 
independent of a2. A general proof valid for arbitrary potential and all cor- 



relation functions was given by Akemann et al. c^. The essence of the proof 
is a remarkable generalization of the identity for the Laguerre polynomials, 
lim„^oo Ln{x/n) — Jo(2y/x) , to orthogonal polynomials determined by an ar- 
bitrary potential. It was proved by taking the continuum limit of the recursion 
relation for orthogonal polynomials. 
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In the second class, an arbritrary fixed matrix is added to W in the Dirac 



operator (107). It has been shown that the microscopicj-snejCtral density and 
the microscopic spectral correlations remain unaffected L30e3 for parameter 
values that completely modify the average spectral density. 

Microscopic universality for deformations that affect the macroscopic spec- 
tral density implies the existence of a scale beyond which universaliti£_breaks 
down. It can be interpreted naturally in terms of the spreading widthEa which 
is the equivalent of the Thouless energy. 

Based on the general form of of the pqChPT partition function one could 
argue that universality of the microscopic correlators in chRMT is automatic. 
However, one really has to show the stability of the effective partition func- 
tion with respect to variations of the distribution of matrix elements. For 
the Wigner-Dysao-ensembles the stability of the saddle-point manifold was 



demonstrated intll 



For reasons of mathematical simplicity, most universality proofs have been 
performed for the (3 — 2 ensembles. Recently, we have attempted to prove 
universality by establishing relations between the kernels for the (3=1 and 
(3 = 4 correlation functions and the kernel for the (3 = 2 correlation functions. 
For the Gaussian ensembles such relations are exact identities. iHl. However, 
for an arbitrary potential they are only valid asymptoticallylHSlHj. 

There are many other results related to the uni-jjjprsality of chRMT. We 
mention recent results in QCD in three dimensions [^If'^^I^P^I, results for Dirac 



operators satisfying the Ginsparg- Wilson relation Ii23, relations between the 
.microscopic spectral density and partition functions with two additional flavors 
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6 Chiral RMT as an Exact Theory for the Fluctuations of Dirac 
Eigenvalues 

6.1 Statistical Analysis of Spectra 

Spectra for a wide range of complex quantum systems have been studied both 
experimentally and numerically (a excellent recent review was given by Guhr, 
Miiller-Groeling and Weidenmiiller 0) . One basic observation is that the scale 
of variations of the average spectral density and the scale of the spectral fluc- 
tuations separate. This allows us to unfold the spectrum, i.e. we rescale the 
spectrum in units of the local average level spacing. Speciflcally, the unfolded 
spectrum is given by 

A^- r\p(.\'))d\', (117) 

J — CXD 
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with unfolded spectral density 

Pu„f(A) = ^<5(A-ArO- (118) 

k 

The fluctuations of the unfolded spectrum can be measured by suitable 
statistics. We will consider the nearest neighbor spacing distribution, P{S), 
and moments of the number of levels in an interval containing n levels on 
average. In particular, we will consider the number variance, T,2{n), and the 
first two cumulants, 71 (?t.) and 72 (?t-). Anothfir useful statistic is the A3(n)- 



statistic introduced by Dyson and MehtalHj. It is related to T,2{n) via a 
smoothening kernel. The advantage of this statistic is that its fluctuations as 
a function of n are greatly reduced. Both S2('t-) and A3(n) can be obtained 
from the pair correlation function. 

Analytical results for all sneiitral correlation functions have been derived 
for each of the three ensembleslHj via the orthogonal polynomial method. We 
only quote the most important results. The nearest neighbor spacing distribu- 
tion, which is known exactly in terms of a power series, is well approximated 
by 

P{S)r^ S/^expi-apS"^), (119) 

where ap is a constant of order one. The asymptotic behavior of T,2{n) and 
A3(n) is given by 

S2(n) - (2/71^/3) log n and A3(n) - /3S2(n)/2. (120) 

Characteristic features of random matrix correlations are level repulsion at 
short distances and a strong suppression of fluctuations at large distances. 
For uncorrelated eigenvalues the level repulsion is absent and one flnds 

P(S')=exp(-S'), (121) 

and 

I]2(n) = n and A3(n) = n/15. (122) 

6.2 Results for Spectral Correlation Functions 

Recently, lattice QCD Dirac spectra were calculated and analyzed by a number 
of different groups. It was found that spectral correlations below a given scale 
are given by chRMT, whereas above this scale, deviations from chRMT due 
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Figure 2: The number variance, S2(n) and the first two cumulants, 71 (n) and 72 (n) as a 
function of n. 
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The first cmnplete Dirac spectra on relatively large lattices were obtained 
by Kalkreuter t2j who calculated all eigenvalues of the N,. = 2 Dirac operator 
both for Kogut-Susskind (KS) fermions and Wilson fermions for lattices as 
large as 12''. In both cases, the Diraojuatrix was tri-diagonalized by CuUum's 
and Willoughby's Lanczos procedureE£D and diagonalized with a standard QL 
algorithm. This improved algorithm makes it possible to obtain all eigenvalues. 
This allows us to test the accuracy of the eigenvalues by means of sum-rules 
for the sum of the squares of the eigenvalues of the lattice Dirac operator (see 
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eqs. ( pH pT])). Typically, the numerical error in the sum rule is of order 10 *. 
Results for the spectral correlations for eigenvalues calculated by Kalkreuter 



!£j both for KS and Wilson fermions are shown in Fig. |[ The simulations 
for KS fermions were performed for 4 dynamical flavors with ma = 0.05 on a 
12^ lattice (the lattice spacing is denoted by a). The simulations for Wilson 
fermions were done for two dynamical flavors on an 8'^ x 12 lattice. 

The results for T,2{n), 71 (n) and ^2{n) obtained by spectral averaging show 
an impressive agreement with the RMT predictions. In the context of RMT 
it has been shown that spectral averages and ensemble averages coincide ll^. 
This property is known as spectral ergodicity. 

Spectra for different values of /3 have been analyzed as well. It is probably 
no surprise that random matrix correlations are found at stronger couplings. 
What is surprising, however, is that even in the weak-coupling domain (/3 — 
2.8) the eigenvalue correlations arc in complete agreement with Random Matrix 
Theory. 

In the case of three or more colors with fundamental fermions, both the 
Wilson and Kogut-Susskind Dirac operator do not possess any anti-unitary 
symmetries. Therefore, our conjecture is that in this case the spectral corre- 
lations in the bulk of the spectrum of both types of fermions can be described 
by the GUE. This was recently confirmed for a wid£-|pjige of /3- values both 
below and above the deconfinement phase transitioncJCj. 

As far as correlations in the bulk of the spectrum are concerned, in the 
case of two fundamental colors the continuum theory and Wilson fermions are 
in the same universality class. It is an interesting question of how spectral 
correlations of KS fermions evolve in the approach to the continuum limit. 
Certainly, the Kramers degeneracy of the eigenvalues remains. However, since 
Kogut-Susskind fermions represent 4 degenerate flavors in the continuum limit, 
the Dirac eigenvalues should obtain an additional two-fold degeneracy. We are 
looking forward to more work in this direction. 



6.3 Correlations near Zero Virtuality 

Spectral ergodicity cannot be exploited in the study of the microscopic spectral 
density and, in order to gather sufficient statistics, a large number of indepen- 
dent spectra is required. One way to proceed is to generate instanton-liquid 
configurations which can be obtained mujch more cheaply than lattice QCD 
configurations. Results of such analysis lla show that for Nc = 2 with funda- 
mental fermions the microscopic spectral density is given by the chGOE. For 
A^c = 3 it is given by the chGUE. One could argue that instanton-liquid config- 
urations can be viewed as smoothened lattice QCD configurations. Roughening 
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Figure 3: The distribution of the smallest eigenvalue (left) and the microscopic spectral 
density (right) for two colors and /3 = 2.0. 



such configurations will only improve the agreement with Random Matrix The- 
puy. Recently, these expectations were confirmed by lattice QCD simulations 
Eil. Results for 1416 quenched SU{2) Kogut-Susskind Dirac spectra on a 10** 
lattice are shown in Fig. ^. We show both the distribution of the smallest 
eigenvalue (left) and the microscopic spectral density (right). The results £2; 
for the chGSE are represented by the dashed curves. 

Agreement of the microscopic spectral density with chRMT for Nc = 2> 
was first demonstEated by means of the valence quark mass dependence of the 
chiral condensate E3. Recentlv—these results were confirmed by the calculation 
of complete Dirac spectra clrcScS. Other recent interesting results are the 
calculation of the micmscopic spectral density for the Dirac operator in the 
adjoint representationCJ and results for the distribution of the smallest Dirac 
eigenvaluCj-at nonzero topological charge obtained by means of the overlap 
formalism O. 



6.4 The Valence Quark Mass Dependence of the Chiral Condensate 

The microscopic spectral density near zero virtuality was first studied in terms 
of the valence quark mass dependence of the chiral condensate. The lattice 
data for S(?7i„) were obtained by the Columbia group l3 for two dynamical 
flavors with sea-quark mass ma = 0.01 and A^c = 3 on a 16'^ x 4 lattice. In 
Fig. 4 we plot the ratio T,{my)/I] as a function of a; = m^VY, (the 'volume' 
V is equal to the total number of Dirac eigenvalues) , and compare the results 
with the universal curves obtained from pq-QCD and chRMT (see eq. (p7|)). 
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Figure 4: The valence quark mass dependence of the chiral condensate S(m) plotted as 
S(m„)/E versus ra^YYl. The dots and squares represent lattice results by the Columbia 
group L3 for values of fi as indicated in the label of the figure. 



We observe that the lattice data for different values of /3 fall on a single curve. 
Moreover, in the mesoscopic range this curve coincides with the random matrix 
prediction for iV/ — v = Q. Of course this is no surprise. In the region where 
valence quark mass is much less than the current quark mass, the fermion 
determinant has no bearing on the Dirac spectrum and we have effectively 
N f — 0. On a lattice the zero mode states and a much larger number of 
nonzero modes states are completely mixed, and therefore it is not surprising 
that the effective topological charge that reproduces the lattice results is equal 
to zero. 

7 Chiral Random Matrix Theory as a Schematic Model for the 
Dirac Spectrum 

7. 1 Introduction 

Up to now we have mainly focusscd on chiral Random Matrix Theory as an 
exact theory for correlations of QCD Dirac spectra on the microscopic scale. 
However, in the literature one finds another important application of RMT, 
namely as a schematic model for a disordered system. Two welliknown ex- 
amples in this category are the Anderson model for localizationlllJ and RMT 
applied to the theory of random surfaces in connection with quantum gravity 
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In this section we will study chiral random matrix theory as a schematic 
model of the chiral phase transition at nonzero temperature and chemical po- 
tential. 

7.2 chRMT at Nonzero Temperature 

In order to obtain a chiral random matrix model for QCD at nonzero temper- 
ature we first separate the Dirac operator in two pieces 

D = D + jodo- (123) 

In a chiral basis with basis functions given by 

^knix) ^ M£)e~'^^"+-'^''° (124) 

the matrix corresponding to 7o<9o is block diagonal with nonzero matrix ele- 
ments given by the Matsubara frequencies. In our model we only keep the 
lowest Matsubara frequencies zLttT. We expect that this is a reasonable ap- 
proximation in the neighborhood of the critical temperature and beyond. The 
matrix corresponding to D is replaced by a chiral random-oiatrix. 
The chRMT Dirac operator at T y^ is thus given byOH 

--'^-(.;.':>(i"^). (-) 

where VIt is given by 

and 1 is a unit matrix. 

If we note that the probability distribution of the random matrix ensemble 
is invariant under W — > UWV~^ and use the equivalent of i<Jicr^(T2 = 1, the 
random matrix Dirac operator can be simplified to 



Notice that the dimension of the identity matrix in eqs. (U2q) and (127) differs 
by a factor 2. 

Our finite temperature chiral random matrix model is thus given by 

Z(z., Nf) = f dWdet'^f {D{T) + m)e-"^"^^^'^^' . (128) 

In the next we will perform a saddle point analysis of this model. 
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1.3 Analytical Solution 



In this section we evaluate the partition function (106) using methods which 
are standard in the supersymnietric formulation of random matrix theoryliililB. 
For simplicity we consider the model for [3 = 2 and choose a diagonal mass 
matrix with equal quark masses. The first step is to perform the average over 
W by performing a gaussian integral. This leads to a four-fermion interaction. 
After averaging over the matrix elements of the Dirac operator, the partition 
function becomes 

Z{i^,Nf) = J Vi;*V^/j cxp[ - —^i}il^R,^R*,i^ik+'n^i'R\^Ri+^^L*k^Lk 

-^7^^(V'^*^{,+^{*V^,)]. (129) 

The four-fermion terms can be written as the difference of two squares. Each 
square can be linearized by 

/2 
dcrexp(-^ -iQcr) . (130) 

The different Q variables, can be combined into a single complex Nf x Nf 
matrix, A, resulting in 

Z{v,Nf)^ f VAV^V^* exp[-^^^:^-^TTAA'' 

- i^il'^lM + ™) - ii^R*i4^Ri{A'^ + m) - i7rT(V'^*V'{ . + V{fcV'L)] • 

(131) 

Note that the temperature-dependent term can be rewritten as 



Using this, the fermionic integrals can be performed, and the partition function 
with |i/| more left-handed modes than right-handed modes is given by 

Z(^, Nf) = JvAeM-^TrAA^]de&\iA + m)det" (^^+ J^ ^^""^^^ 

(133) 

Here, A is an arbitrary complex matrix and m is proportional to the identity 
matrix. 
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The chiral condensate is defined by 

m=-^d^\ogZ , (134) 

For n — > cxD it can be evaluated with the aid of a saddle point approximation. 
The saddle point equations are given by 



n/3S^ 



A + n{A + m){{A'' +m){A + m)+TT^T'^) ^=0 . (135) 



An arbitrary complex matrix can be diagonalized by performing the decompo- 
sition 

A^UAV-\ (136) 

with all eigenvalues positive and U and V unitary matrices. We find that the 
solution of (135) yields U ^ V ~ 1 with eigenvalues A^ given by the positive 
roots of 

T.^Ak{{Ak + mf+7T^T^)-Ak-m = . (137) 

In order to calculate the condensate, we express the derivative of the partition 
function in (3.10) in terms of an average over A, 

\^^^^\-2^^^^{-t^T'A^y^ • (i3«) 

Below Tc and for to — > we find from the saddle point equation, 

|(l/;V')| -S(l-7r2r2£2-)l/2 ^ ^^gg^ 

In the chiral limit we thus find a critical point at 

Tc^^ ■ (140) 

Note that the dimensions are matched after including the mode density N/V = 
1. Here, and elsewhere, our convention is that N = 2n. At Tc the solution of 
the saddle point equation develops a non-analytic dependence on m resulting 
in the condensate 

|(-0V)| = S]3mi . (141) 

Therefore, we reproduce the mean field value for the critical exponent <5 = 3. 
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It is also possible to obtain an analytical expression for the spectral density 
of this model as a function of the critical temperature. It can be obtained by 
solving the Schwinger-Dyson equations for the resolvent of the modelEll or from 
the observation that the solution of (135) does not depend on the number of 
flavors. For Nf -^ the spectral density is then obtained from the discontinuity 
of the chiral condensate across the real m-axisEj. The spectral density thus 
follows from the solution of a cubic equation. The zero temperature limit is 
a semicircle whereas the high temperature limit is given by two semi-circles 
separated by 2ttT. 



7.4 chRMT at Nonzero Chemical Potential 

The chemical potential enters in the continuum QCD partition function as 



170^0 -^ ^70^0 + fJ-Jo- 

resulting in the chRMT Dirac operator E3 

iW + i-i 




^(A')=UM/t+^ 



(142) 



(143) 



If we nnfipp that fijQ has the same reality properties as 7o9o, we immediately 
concludell£j that these ensembles can be classified according to the same anti- 
unitary symmetries as the chRMT's for fj, — 0. However, the term proportional 
to fj, violates the anti-Hermiticity of the Dirac operator, and typically the 
eigenvalues of D(fj,) will be scattered in the complex plane. 

Many of the problems associated with the presence of a chemical potential 
are related to the loss of non-Hermiticity. Exactly, these problems are repro- 
duced by the chRMT partition function with Dirac operator given by (143). 
For examji 
imationEf^ 



^ ir this model one can study the problems pi-tl-i^ quenched approx- 
I43|ji44| ^YiQ structure of the Yang-Lee zeros EfjEll and the problems 
with the Glasgow methodlil^. Below we will discuss the first two applications. 



7.5 Yang-Lee Zeros 

The QCD partition function 

Z(to, /i) = / det^^ {D + m + fijo)e~^^'" 



(144) 



is a pnlynnmial in m and /i. For simplicity, let us consider the case m — 0. 
Then 
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ZO^^(m = 0,Ai)^n(^-^fc) 



(145) 
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and the baryon density is given by 

nB = ^d^,\ogZ{m^0,l^) = ^J2^^- (146) 

k 

With the zeros scattered in the complex /i-plane we can interpret the real and 
imaginary parts of n^ as the electric field at /i due to charges located at ^k ■ 

Physically, we expect that hb = for /i < ^^ t^ and jumps to a finite 
value at fj, — lie- We thus expect a first order phase transition at fXc- Two 
possibilities for the distribution of the zeros come to mind. First, according 
to Gauss law, a homogenous distribution of zeros fXk along the complex circle 
|/i| = lie results in a zero baryon number density for |^| < ^c- (Of course in the 
thermodynamic limit, a subleading number of zeros may be present inside the 
circle.). Second, we all know that the electric field is zero between two infinite 
parallel plates with equal constant charge density is zero. The analogue of 
this phenomenon in two dimensions is that the electric field is zero between 
two parallel infinite line charges with constant charge density. Therefore, the 
second possibility is that the zeros in the complex ^-plane are located at 

fik = ±l^c + ^, kodd. (147) 

Zp 

We also expect that the baryon number density increases smoothly for fi > ii,,- 
This requires a radially symmetric cloud of zeros outside of the radius |/i| — ij-c 
in the first case, and a constant density parallel to the imaginary axis for 
|Re(^)| > He in the second case. 

The simplest possible model for spherically symmetric distributed zeros of 
the partition function, is one with zeros given by the roots of unity, 

A*fe-Mce^, k = l,---,N. (148) 

This results in the partition function 

k 

For the baryon number density we then obtain 



1 a i.„^_ M 



Af-l 



ns = -rrdalogZ 



i for |m|>Mc fo,^^^. (150) 

for |/i| < He 
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In the second case with zeros given by (147), the partition function is given 
by 

^(a^) = n (^ - ^- - ^)(^ + ^- - 5)- (151) 

k odd 

Up to a constant the zeros are given by the zeros of the cosh-function. We 
thus find that 

Z{p) ^ cosh(/3(^e - ^)) cosh(/3(/ie + A^)) 
ZiO) cosh2(Mc) ' 

The baryon density is then given by (in an arbitrary normalization) 

J_ z( \= sinh(2/?/i) 

2^0-^ og ^^j cosh(2/3/i) + cosh(2/3^e) 

For /3 — > cx) this can be approximated by 



UB = — O^logZ(^) = _^,^^„, , l„Lo«.. ^ ^ (153) 



^B = ^^^2m,.-^y (154) 



and we find that 



ub = for /i < /ic, (155) 

ub = 1 for /i > ^c- (156) 

In the next section we wih study the baryon number density in the chiral 
random matrix model at nonzero chemical potential. 

7.6 Phases in chRMT at ^ ^ 

The chRMT partition function at /i 7^ is obtained from the finite T partition 
function by the replacement T — > i/i. For Nf = 1 this results in the a model 

Zi^i)= j dcjda*{\a\''-p?fe~''\"\\ (157) 



For iV — > 00, the integrals in this partition partition function can be performed 
by a saddle point approximation. We find that 

(7 = for ^i> ^ic ^ Z = ^2^, (158) 



CT = ^1 + ^2 foj. ^<.^^ ^ Z = e-^^^''+^\ (159) 

The critical point is given by the trancendental equations /i^ = exp(— 1 — /i^) 
which is solved by /Ltc ~ 0.53 
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Figure 5: The zeros of the partition function in the complex /i plane (left) and the result 
obtained from a mean field analysis (right). Results are given for for m = (upper), 
m = 0.30 (middle), and m = 1.0 (lower) for N = 192. The zeros of the discriminant of the 
cubic equation are denoted by stars. Note that the scale on the x-axis of the lower figure is 
different. 



We find that the baryon number density above the critical point shows the 
same behavior as in our naive model. However, the baryon density for fj, < fic 
is nonvanishing which disagrees with our expectation for QCD at zero tem- 
perature and finite density. The remedy should be clear. In order to obtain a 
sha rp Fermi-Dirac distribution one has to sum over all Matsubara frequencies 
1113. This point was ignored in our model. In the next section we discuss a 
RMT CT-model that explicitly includes the lowest two Matsubara frequencies. 
This model can be trivially extended to include all Matsubara frequencies, and 
the infinite product jCan be evaluated analytically in terms of a cosh- function. 
For /3 — > oo we findE3 that the zeros of the partition function in the complex 
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/i-plane are located on the line Re(/i) = ±fic- 

A much more disturbing observation is that in the thermDjdynamic limit 
the partition function for ^ < Hc becomes exponentially smallElfEy. This expo- 
nential cancellation is due to the phase of the eigenvalues. It makes accurate 
numerical simulations forbiddingly difficult. 

The random matrix partition function is a polynomial in m and ^ and it 
is straightforward to study the zeros of the partition function in the complex 
TO-plane or the complex /i-plane. To assure our numerical accuracy we have 
determined the zeros of the polynomials by means of multi-precision algorithms 

using as much as 1000 significant digits. The results are shown in Fig. 
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5. The branch points at the end of a cut are the points where two saddle- 
point solutions coincide. In the present case they are given by zeros of the 
discriminant of a cubic equation. The line of zeros can also be obtained from a 



saddle point ana. 



vais (see right half of the figure). For more details and results 



in other fields t^'E^ll we refer to the original literature. 



7. 7 Failure of the Quenched Approximation 

For QCD with three or more colors and fermions in the fundamental repre- 
sentation, the fcrmion determinant is complex. The presence of this phase 
makes Monte-Carlo simulations impossible. A way out might be the quenched 
approximation, i.e. ignore the fermion deterrninant. altogether and hope that 



it will work. This approach was followed intlStlj with the conclusion that 
the critical chemical potential is given by the pion mass instead of the nucleon 
mass. Obviously, this result is physically wrong! 

This problem was first analyzed in chRMT by StephanovEil. He showed 
that the quenched limit is the limit Nf ^ of a partition functions with 
fermion determinant 

|det(i:>(/i)-f to)|^^ (160) 

rather than 

{det{D{^x)+m))^f. (161) 

The absolute value can be represented as (let us take Nf = 2 for simplicity) 

det(7^(/z) + m) dei{D\y) + vi*) = j D^D^'e^^°^f''>+"''>^+^^^°'^''^+"''^^\ 

(162) 
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and Goldstone bosons made out of quarks and conjugate anti-quarks carry a 
net baryon number resulting in a critical chemical potential given by the pion 
mass. 



The random matrix model corresponding to (162) can be mapped onto a 
nonlinear cr-model which can be solved by a saddle-point approximation. The 
exact solution of this model confirms that the critical chemical potential is 



proportional to \pn, the same dependence as for the pion mass. 



1.8 The Chiral Phase Transition in chRMT 

In previous section we have studied a model at nonzero temperature and a 
model at nonzero chemical potential. These two models can be combined in 
the following schematic chRMT model E3 for the chiral phase transition, 

Z= f DWdct^^ ( ,„, ™ .^ '^ + '^) (163) 

with C a diagonal matrix with Cu — airT — ibfi for one half of the diagonal 
elements and Ck ~ ^anT ~ ib^ for the other half with a and b dimension- 
less parameters. This model can be considered as the matrix equivalent of a 
Landau-Ginzburg functional. We thus expect to find mean field critical ex- 
ponents. The advantage over using Landau-Ginzburg theory is that in this 
case the spectrum of the Dirac operator is accessible. In particular, this might 
reveal interesting behavior of the Dirac spectrum near the critical point. For 
example, for /i = we have shown that the fluctuations of the smallest eigen- 
value of the Dirac operator can be used as an order parameter for the chiral 
phase transitional. 

Also in this case the partition function can be reduced to a cr-model. For 
the simplest case of Nf = 1 it is given by 

Z{T,fi)^ f dae-^^^''\ (164) 

where 

n{a) = aa"^ - log((cr + m){a'^ + m) - {fi + iTf) 

~ log((a + TO)(o-t + to) - (/i - iTf). (165) 

The chiral condensate is given by the expectation value of a. The saddle-point 
equation is a fifth order equation in g. Therefore this model is very similar to 
a i/)^ Landau-Ginzburg theory. For tti = we find that a is real and the saddle 
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point-equation is given by 

a[(T^ - 2(^2 _ y2 ^ 1^2 ^ (^2 ^ y2)2 ^ ^2 _ y2] ^ ^ (^gg) 

The critical points occur where one of the solutions of the cubic equation 
merge with the solution cr = 0, i.e. along the curve (/i^ + T^)^ + ^^ — T^ = 0. 
At the tri-critical point three solution merge. This happens if in addition 

In Fig. 6 we show the phase diagram in the fiTm space. In the m = 
plane we observe a line of second order phase transitions and a line of first 
order phase transitions. They join at the tricritical point. Also joining at the 
tricritical point is a line of second order phase transitions in the 77i-directions 
which is the boundary of the plane of first order transitions in yLtT?Ti-space. 




0-4 f^ 0.5 



Figure 6: Phase diagram of QCD with two light flavors of mass m as calculated from 
the random matrix model. The almost parallel curves on the wing surface are cross 
sections of this surface with m =const planes. 



7.9 Possibility of a Localization Transition in QCD 

From the theory of Anderson localization we know that the eigenvalues corre- 
sponding to localized states are statistically independent and thus obey Poisson 
statistics. Since it has been well established that in QCD at zero temperature 
the correlations of the eigenvalues of the QCD Dirac operator are given by 
chiral Random Matrix theory, we conclude that the eigenstates are extended. 
This also follows from a direct analysis of the Dirac wave functions for gauge 
field configurations given by a liquid of instantons. However, we wish to point 



out that (Studies with Wilson Dirac fcrmions indicate that cigcnfunctions are 



localizedE£il. We do not have an explanation for this discrepancy. 

One reason to expect a localization transition with increasing temperature 
is dimensional reduction. States are more likely to be localized in lower di- 
mensions. For example, in a one-dimensional disordered system all states are 
localized. One the other hand, because of asymptotic freedom the coupling 
becomes weaker at higher temperatures. Let us analyze more closely what will 
happen. 

Lattice simulations for the valence quark mass dependence of the chiral 
condensate weje performed for temperatures both below and above the critical 
temperature 113. It was found that, in the ergodic domain, all data can be 
rescaled onto a universal curve given by chRMT. This shows that the states 
are extended below the critical point. 



According to a theoretical argument due to Parisit^ given in the context 
of disordered systems a localization transition can only occur in quenched 
theories. In QCD, this argument can be translated as follows. Let us assume 
that the eigenfunctions of the Dirac operator are localized. Then the joint 
eigenvalue distribution factorizes in one-particle distributions (A^-t-?Ti^)^-''F(A). 
The average spectral density is thus given by 

/n 
d\2--- dK \{{\l + m^f!F{X^) ■ ■ ■ F{K) 
fe=i 

= ci(A2 + m2)^/i^(Ai). (167) 



We then find 



lim lim lim p{\i) = 0, (168) 



i.e., no spontaneous breaking of chiral symmetry. We thus conclude that in the 
chiral limit no localization can occur in QCD with light quarks. What happens 
in quenched theories remains an open question. However, studies with long 
range interactions given by ^ ±l/|i?i — Rj\'^ with random signs and random 
positions Ri as is the_case for instanton liquid simulations indicate that all 



states are delocalizedli^ 



7.10 Triality at ^i ^ Q 

Let us finally discuss how Dirac spectra at nonzero chemical potential depend 
on the Dyson index of the matrix elements. We remind the reader that the 
classification of Dirac operators with chemical potential is the same as for zero 
chemical potential. 
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Numerical simulations have been performed for all three classes. A cut 



•y ayi« below a cloud of eigenvalues was found in instanton 
1£3 for Nc = 2 at fj, y^ which corresponds_±D /? = 1. In 
lattice QCD simulations with staggered fermions for Nc — 2 



along the imagina: 
liquid simulations 



' P^^'^ a depletion 
of eigenvalues along the imaginary axis was observed, whereasJior Nc = 3 the 



eigenvalue distribution did not show any pronounced features 
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Figure 7: Scatter plot of the real (x), and the imaginary parts (y) of the eigenvalues of 
the random matrix Dirac operator at nonzero chemical potential. The values of /3 and fj, are 
given in the labels of the figure. The full curve shows the analytical result for the boundary. 



In the quenched approximation, the spectral properties of the random 



matrix Dirac operator (143) can easily be studied numerically by diagonalizing 



a-Sfit of matrices with probability distribution (FTOg) . In Fig. 7 we show results 
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for the eigenvalues of a few 100 x 100 matrices for fi = 0.15 (dots). The 
solid curve represents the analytical result for the boundary of the domain of 
eigenvalues dedved inO for (3 = 2. However, the method that was used can 
be extendedll£3 to (3 = 1 and f3 = A and with the proper scale factors we find 
exactly the same solution. 

For (3 = 1 and (3 = 4 we observe exactly the same structure as in the 
previously mentioned (quenched) QCD simulations. We find an accumulation 
of eigenvalues on the imaginary axis for (3 = 1 and a depletion of eigenvalues 
along this axis for (3 = A. This depletion can be understood as follows. For 
fi = all eigenvalues are doubly degenerate. This degeneracy is broken at 
/i ^ which produces the observed repulsion between the eigenvalues. 

The number of purely imaginary eigenvalues for /3 = 1 scales as \/N and 
is thus not visible in a leading order saddle point analysis. Such a yN scaling 
is tjipjjcal for the regime of weak non-hermiticity first identified by Fyodorov et 
al. 
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Using the supersymmetric method of random matrix theory the yN 
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dependence was obtained_analytically by EfetovE£3. Also the case /3 — 4 was 
analyzed analytically inE£3 with results that are in complete agreement with 
our numerical simulations. Obviously, more work has to ha, done in order to 



arrive at a complete characterization of universal featuresll£3 in the spectrum 
of nonhermitean matrices. 

8 Closing Remarks 

In these lectures we have presented analytical results for the infrared limit of 
the QCD Dirac spectrum. We have shown that the correlations of the Dirac 
eigenvalues for level spacings below the Thouless energy are given by the zero 
momentum limit of the partially quenched chiral partition function and agree 
with results obtained from chiral random matrix theory. A key ingredient in the 
formulation of the effective theory is the structure of the integration manifold. 
We have argued that it is given by a super-Riemannian manifold which is 
characterized by a symbiosis between compact and noncompact degrees of 
freedom. We have shown that the valence quark mass dependence of the 
chiral condensate that follows from the zero momentum sector of the effective 
chiral partition function coincides with the result from chiral Random Matrix 
Theory. In this formulation, universality is natural. The structure of the chiral 
Lagrangian only depends on the pattern of chiral symmetry breaking. 

An important condition for the validity of the chiral Lagrangian is the 
requirement that the only low-lying modes are the Goldstone modes associated 
with the spontaneous breaking of chiral symmetry. Without confinement this 
would not be the case. On the other hand, according to the Bohigas conjecture, 
eigenvalue correlations are given by random matrix theory if the corresponding 
classical system is chaotic. If the classical motion of quarks in 4+1 dimensions 
is chaotic the eigenvalues of the Dirac operator are correlated according to 
chRMT. Then necessarily the low-energy dynamics is given is given by the 
chiral Lagrangian and this can only happen if we have confinement. 

Can we reverse this statement? Is the classical motion of a fermion in a 
confining theory necessarily chaotic? The point I wish to make is that the 
Bohigas conjecture and confinement are not unrelated. For a deeper under- 
standing of confinement a proof of the Bohigas rnj^jpr|:nr^ ^^z\^ be required. 
Recent progress in this direction is encouragingE 
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